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strands of hydrogen bonds constituting the a— helical backbone. 
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1. Introduction 

Flexibility of proteins is the feature giving rise to some of their remarkable properties 
[1]. Even minute conformational changes can modify long-range electronic interactions 
in proteins and open pathways for molecular motions, provided proteins are not rigid 
object |2]. In this context, we raise the following questions: Is the flexibility of proteins 
relevant for energy storage and transfer? Can local protein conformation come in 
discrete portions of elastic energy? Is this mechanism quantum or not? The ubiquitous 
a— helical sections of the secondary structure of proteins have already been proposed 
by Davydov and co-workers as a conduit for quantum coherent energy transfer due to 
the interchain hydrogen bonded network which stabilizes their structure and forms a 
path along the a— helix [3]. Here we address these questions from the stand point of 
measuring the conformation of the a— helix through a set of natural geometric quantities 
such as curvature and torsion and their coupling to a quantum of vibrational energy 
present on the helix. To answer the above questions, we propose a nonlinear quantum 
mechanism of energy transport due to a geometry coupling effect along the hydrogen 
bonded lattice of a— helical sections of proteins. But first let us introduce the basic 
structure of proteins. 

H H 
I I 

HO— C— C— N— H 

II I 
R 

Figure 1. Chemical structure of an amino acid, where the symbol R denotes the 
residue, a side hydrocarbon chain. 

The primary structure of proteins is a chain sequence of amino acids, biomolecules 
which have the structure depicted in Figure [TJ Amino acids are bound together, to form 
proteins, by a peptide bond constituted by the elimination of a water molecule from 
the acidic group and the — NH terminus. To each amino acid is attached a residue, 
a side hydrocarbon chain. There exist 20 different residues and their sequence gives 
the protein's primary structure. Apart from these residues, proteins are constituted 
by a repeated pattern (see Figure |2|: a group of atoms, called the peptide group (see 
Figure |2]). The peptide group has two functions: i.) act as a bond between carbon 
atoms holding the distinct residues; ii.) due to their dipole moment make the protein 
fold into secondary structures stabilized by the formation of hydrogen bonds between 
the peptide groups. The two most common distinct shapes found in protein secondary 
structures are a— helices, shown in Figure [sj and planar structures called /3— sheets jl]. 
The mechanism of energy transfer and localization in proteins proposed here concerns 
the a— helical sections. 

The underlying mechanism describing the local conformational dynamics of proteins 
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Figure 2. A chain of amino acids, linked by peptide bonds (peptide groups are 
given in blue boxes), with different residues Ri and R2. Here the peptide groups are 
interconnected by a single C atom holding the residue. The bonds C-N and C-C can 
rotate, leading to a rearrangement of the chain in space. 



for energy transfer motivated the idea of the conformon, a concept advanced by various 
authors [SI El El E] to play an important role in biology. It was originally envisioned 
as an explanation of how the energy 0.5 eV) released in the hydrolysis of the ATP 
molecule, an amount insufficient to excite the electronic structure, is transported along 
biomolecules without being immediately dissipated to the environment's degrees of 
freedom. This energy is slightly larger than two quanta of the Amid-I vibration of 
the peptide group (Amid-I is the C=0 stretch, which can drive the hydrogen bond 
between two peptide groups). The Amid-I vibrational quanta are thus accepted as 
the standard initial sinks in the biophysical processes utilizing the energy released in 
the hydrolysis of the ATP molecule [3]. The conformon is a concept spanning a bridge 
between the three physical processes in the base of the functioning of biological molecules 
such as proteins: i.) energy storage and transfer; ii.) conformational changes; iii.) 
charge transport. The conformon was introduced independently by Green and Ji [5], by 
Vol'kenstein [B] and by Kemeny and Goklany [7] . Green and Ji proposed the conformon 
as "the free energy associated with local conformational strain." Vol'kenstein defined 
it as "the displacement of the electronic density in a macromolecule that produces the 
deformation of the lattice, i.e. the conformational strain." Kemeny and Goklany, using 
Holstein's theory [10], defined the conformon as an electron spread over a molecule 
forming a polaron bound to a number of phonons. This bound state can then propagate 
from molecule to molecule. 

To avoid confusion it should be mentioned that the concept of conformons as 
bosons of linear relaxation modes with photon statistics is also used in the context 
of polymer networks and dissipative liquid-like system such as glass forming fluids [TT] . 
It is interesting to note that the conformon in the definitions presented above is a 
quantum phenomenon and is subject to the field of quantum biophysics recently covered 
in [12]. Other conformon-like concepts in biophysics have been introduced [131 [E] but 
are classical in nature. 

Here we propose a physical definition of the conformon as a nonlinear mechanism 
of exchange of energy and coherence between the quantum of vibrational excitation and 
the conformational degrees of freedom of the a— helix which, due to the presence of the 
excitation, are quantized themselves. By introducing a measure for the conformation 
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of the a— helix, we will show that the conformon can be viewed as a quantum of local 
conformation. 




Figure 3. Representation of tlie (3D) lielical structure of tire a— helix portion 
of a protein. The peptide groups are represented by circles. Hydrogen bonds 
(represented by the arrows) link every third peptide group. This forms the three 
strands (represented in three different colours) constituting the backbone of the helix. 

We start the description of the conformational dynamics of biopolymers with their 
representation as thin rods where the relevant degrees of freedom are the torsion and 
the curvature. In the Kirchhoff model of an inextensible rod with fixed torsion, the 
curvature dynamics allows for solitary waves in the case of non-circular cross-sections [9] . 
A quantum particle constrained in a rod is subject to an attractive potential induced by 
the torsion and the curvature [15]. Applying this to the problem of electron transport, 
it can be shown that this quantum potential leads to localized states, as the energy 
density in the polymer is proportional to the curvature squared. This is essentially 
because the rod's curvature and torsion "interact" with the excitation by inducing a 
potential well, which traps the electron, in addition to creating a nonlinearity in its 
Schrodinger equation |9]. While this scenario is somewhat analogous to a polaron (a 
localized electronic bound state in a discrete lattice which is not perfectly periodic) 
the origins of the two mechanisms are quite distinct, with the curved geometry of the 
biopolymer playing a key role in the creation of the conformon. 

Our motivation for introducing the notion of the conformon is its potential for 
an explanation of the recently observed long-lived low frequency nonlinear vibrational 
states in predominantly a— helical proteins. Pump-probe experiments revealed that 
low frequency nonlinear modes are essential for functionally important conformational 
transitions in proteins containing a— helices [23] • The lifetime of these states for 
bacteriorhodopsin can go over 500 ps, which is intriguing considering the expected 
dissipation and dephasing timescales for proteins of that size. Other nonlinear 
mechanisms of energy transport in a— helices, as Davydov's theory of extrinsic 
anharmonicity of hydrogen bonds coupled to an excitation localized on the lattice, fail to 
exhibit such a long lifetime [HI [IS] considering the physical parameters of the hydrogen 
bonds and couplings in the models. 

While the Dandoloff-Balakrishnan model P] just discussed provides a possible 
dynamical underpinning for the conformon, the aim of the present paper is to give 
the underlying microscopic quantum mechanical description. This task is achieved by 
considering the full structure of the hydrogen bonded network of the a— helix (see Figure 
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[3]). In reality, a— helices are fixed by three strands of soft hydrogen bonds stemming 
from the periodically placed peptide groups. Except for reference |T6], previous studies 
simplify the description by assuming an effective single strand of hydrogen bonds, thus 
concentrating on the helical-symmetry preserving cases. Here we address the symmetry 
breaking modes analytically, which not only have been numerically demonstrated to be 
less energetic and thus more stable [H [16], but also give the realistic description of the 
a— helix hydrogen backbone. 

In the following section we construct the hamiltonian for the energy transfer along 
a— helical proteins taking into account the three strands of hydrogen bonds. Next, in 
section [3| we solve the proposed theory. Section |4] discusses the conformon solution and 
justifies it as a quantum of conformation to be found in a— helical proteins. 



2. The Hamiltonian 



The a— helix of a protein contains u = 3.6 peptide groups per period of the helix[3]. 
The equilibrium position of each peptide group is given, in a coordinate system centered 
with the axis of the helix (see Figure [s]), by the radius vector 



f2TTi\ ^ . f2TTi\ 



Cxcos + e„sm , 

\ y J \ ^ /\ ^ 



+ e. — , (1) 



where 6^,6^ and Cz are the orthonormal basis vectors, i is an integer labeling each 
peptide group, r = 1.7 A is the helix's radius and a = 5.4A is the period of the helix 
along z Along the a— helix, a peptide group i is connected to its neighbours i — 1 
and £ + 1 by rigid covalent bonds. The same peptide group i is bound with its £ — 3 and 
i + 3 neighbours through a soft hydrogen bond. These constitute the three hydrogen 
strands along the helix (see Figure |3]). 

Due to the softness of the hydrogen bonds, the peptide groups can be displaced 
from their equilibrium positions. Let us consider sets of three consecutive peptide groups 
(along the helix): according to Euclidean geometry, each set will determine an unique 
triangle/plane containing the respective peptide groups. We can now look at the a— 
helix as a sequence of such triangles/planes connected to each other in two different 
ways: i.) there is a hard covalent bond connecting the last peptide group of one triangle 
to the first peptide group of the next triangle; ii.) each of the three peptide groups of one 
triangle are connected with the peptide groups of the other triangles through the three 
hydrogen bond strands. From the mechanical point of view we may treat each plane n, 
without loss of generality, as a rigid disk with an effective mass and an inertia tensor 
In which take into account the residues attached to each peptide group. In what follows, 
we will model the a— helix as a linear chain of such rigid disks elastically connected with 
three springs representing the hydrogen bonds (see Figure |4]). The covalent bonds will 
be treated as rigid constraints reducing extra degrees of freedom (and thereby are not 
depicted in Figure |4]). To quantify the local conformation in this model of the a— helix 
we must introduce a measure of the relative orientations of the disks. To achieve this we 
track the change in the orientation of the unit normal vector to each disk with respect 
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to its neighboring one. Considering the set of all those vectors, represented one after 
the other and separated by an effective constant distance d (the a— helix is assumed 
inextensible, as discussed later), we construct a discrete space curve F (see Figure [sj. 
Next, we associate elastic energy to every configuration of F, since every deviation from 
the equilibrium position either presses upon or stretches the hydrogen bonds from the 
three strands. This way the conformation of the a— helix in space can be quantified. 




Figure 4. Our model of the a— helix, as a sequence of rigid disks (consisting of three 
consequtive peptide groups and the hydrocarbon residues attached to them) connected 
with three springs represented by the wavy lines corresponding to hydrogen bonds. The 
covalent bonds connecting and fixing the mutual orientation of two consequtive disks 
are treated as contraints. 

We describe the configuration of the a— helix as a discrete curve F comprised of 
nodes (representing disks or 3A^ peptide groups) and N — 1 straight edges of length 
d together with an assignment of local material orthonormal frames in the standard 
Darboux-Ribaucour notation = (&n, '^n, ^n)"^ per node, where T denotes transposition 
[T7] . Here is the unit normal vector to a disk n giving the orientation of the associated 
edge (see Figure |5]). 




disk n-l disk n-t-l 

Figure 5. The discrete space curve F is constructed from straight lines of constant 
length d (called edges) along the normals to the rigid disks making up the a— helix. A 
local orthonormal frame is attached to every node of F. The change of the orientation 
of one such frame with respect to its neighbouring one along F is a measure of the 
elastic energy stored in a particular conformation of F. 

Now, we assign an elastic energy E{T) to any conformation of F. Note, we assume 
the a— helix is inextensible, therefore we do not include stretching energy. It is 
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straightforward to drop this assumption by also including a stretching term. Such a 
term is going to generate an acoustic mode which is not the focus of the present study. 

To compute the elastic energy E{T) we use a three dimensional symmetric]!] strain 
tensor S. In the first approximation, as in Hooke's law, -E(r) is quadratic in the 
displacements from the initial equilibrium configuration as measured by the Darboux 
vector The Darboux vector gives the relation between fa and fa+i- The sequence of 
Darboux vectors /i, /2 • • • /n gives the evolution of V [TH]. Finally, for the elastic energy 
we have: 

1 ^ T 

= 2 E (^"^n - ^"^n ) ^(fn-e), (2) 
n=l 

where E,^ is the initial configuration of F and the Darboux vector in the local material 
frame /n is given by: 

= + + rj^, (3) 

where is the geodesic curvature, k„ is the normal curvature and is the torsion 
defined on each disk of F. In the continuum limit, where N ^ oo and — )■ 0, one 
obtains the bending and twisting energies of the inextensible Kirchhoff rod [19j. 

To compute the curvatures and torsion present in equation ^ we start by 
establishing a mapping between the local material frames /„ and the fixed laboratory 
frame e = (e^, Cy, Cz)^. Then we give the orientation of each frame /n in terms of the 
Euler angles: 

fn = Rn{lpn,On,^n)e, (4) 

where is the product of the three Euler matrices representing three rotations by 0n 
degrees around the 2r— axis, by 9^ degrees around the x— axis and by V'n degrees around 
the 2;— axis again. In this way, the transformation of the frame from disk to disk, i.e. 
the discrete analog of the Darboux- Ribaucour equations [18], are given by 

fn+l = R-n+lRn^ fn, ^ (^fn+1 — fnj = ^nfn, (S) 

where = (-Rn+i-Rn^ — I)/d. In the continuous limit, when N ^ 00 and d ^ 0, 
the elements {Aa)3i, (v4n)23 and (An)i2 of the matrix give, respectively, the geodesic 
curvature, the normal curvature and torsion of the space curve [18]. We compute A^ in 
our discrete case and identifying the same elements we finally obtain the coefficients of 
the Darboux vector in equation (]3|: 

= sm %pa sm fc'n 2 ^ '^'i 2 ' ' 

= COS?/'nSm6'n SlUtpr, , (7) 

Tn = h cose*!! J . (8) 

a d 

I in the the material frame 
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The interpretation of the Euler angles is the following: 6^ and (p^L are the polar 
and azimuthal angles, respectively, of tn in the lab frame e. The angle ip^ rotates the 
(^n, &n) plane around tn and thus leaves tn unchanged. This means that the normal 
vectors are mapped onto the z— axis of the lab frame through precession 0n around 
this axis and tilting 9^ off this axis. The differences in the polar and azimuthal angles 
of t-a between two neighboring disks are not associated with the choice of a cartesian 
system e and are a measure of the bending and twisting of the initial conformation of 
F; ip^ constitutes the gauge freedom in the problem which we fix due to the existence of 
a covalent constraint to each disk. We set ip,^ to be a constant for every n and we also 
fix the strain tensor 

S = a\ a G M (9) 

in the basis /n corresponding to an isotropic elastic response to any change of 
conformation of F. 

The total hamiltonian describing the quantum dynamics of a quantum of vibrational 
excitation present on the a— helix modeled as F is a collection of an elastic part Hy , an 
excitation part H^x and an interaction term Hint between them. A kinetic energy term 
of the form Hkin = \ Xln '^n-^n^^n (wherc 1^ is the moment of inertia tensor of each disk 
and cJn is its angular velocity) can be included but we will omit it in accordance with 
the adiabatic approximation. This is justified due to the enormous mass and moment 
of inertia tensor of the material disks. As a result the excitation's motion adjusts 
instantaneously to the conformation of F, that is for the excitation's dynamics the slow 
dynamics of conformational change is effectively static. 

We are now in a position to write the quantum energy operator for the elastic part 
of the total hamiltonian making the variables from equation ([2]) operator valued. Here 
we choose to define the initial condition for F by setting E,^ = corresponding to the 
equilibrium configuration of F for which the protein is helical. The energy consumed in 
the deviation from this equilibrium state is obtained by substituting ([9]) in ([2]) and then 
using (§, and g: 

n 

n 

The energy operator for the vibrational excitation present in the a— helix Hf,^ is 
given in the standard manner 



Hex = ^ l^-f^O-Bii-Bn — -/(-Bn-Bn+l + B\^Bn-l) 



'111 



where -Bj and are the excitation's creation and annihilation operators for disk n; J is 
the hopping amplitude kept constant due to the inextensibility condition on F. Here, for 
simplicity, we assume the chain is infinite and do not consider the effect of an excitation 
entering the chain from either end or finite boundary conditions. 
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The interaction term Hint between the excitation and orientational degrees of 
freedom is deducible from the inextensibihty condition on F and the local modification 
of the on-disk energy for the curved conformation of F in the nearest neighbor 
approximation, that is E^, ^ Eq + a(g^^f_g^) (^n+i - On) + o(0„^f-^„) (0n+i - 0n)- Finally, 

H^nt = Yl - 0n) ^t^, + Xe{0n+1 " ^n)^!^^} , (12) 

n 

where X4> ^^^^ Xo are coupling constants, small as compared to the conformational energy 
'^d^X(t>/(^ < 1 and 2cPxe/oi < 1. 



The total hamiltonian is the sum of (10), (11) and (12). We explore the theory 



H = Hr + Hex + Hint- (13) 
3. Solving the Theory 



Our analysis of (13) will be based upon a product trial wave function |\E') = |0) in 
which describes a single excitation in F as 

\^) =J2<^n{t)Bi\0)e., (14) 

n 

where |0)ea: is the vacuum state of the its oscillators and |0) is a coherent state for the 
angular displacements, for which the following hold 

{<P\eM=On, (0|4|0)=0n. (15) 

Here a^it) is a complex number representing the probability amplitude for finding 
the excitation in a particular site; 6'n and 0n are the average values of the angular 
displacements of F's edges. Next we will minimize the average value of H with respect 
to the product wave function by calculating H = (\E'| |\E') , where we have to keep in 
mind the normalization condition Xln'^n'^n ~ 1- -^^^ averaged energy function in the 
adiabatic approximation we obtain 

0On«n ~ J{0'n(^n+l + «n'^n-l) + (16) 

n 

+ [Xcf>{4>n+1 - 0n) + Xe{6n+1 " ^n)] + 



a 



m ) 



Next we find the corresponding Hamilton's equations 

ih-^ = {Eo + W + Tn) an - J(an+i + a^-i), (17) 

where Tn = x<^(0n+i — 4>n) + Xe(^n+i — ^n) and the total conformational energy is 
a 

2^2 
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The other two Hamilton's equations for 0n and 6*11 constitute the static (since we are in 
the adiabatic approximation) configuration of F : 



a 



One possible solution is 



2d\e 



a 



'(|an- 



2 

with Tn = —V \ a-a\ and 
a 



2(Px<t> l«n| 



a 



-, — ^n+1 



X2 \an\ 



2<Pxe I On I 



a 



A gauge transformation of the wavefunction (14) 

an(t) = qn{t) exp | -^t (Eq + - 2 J) 



(19) 
(20) 

(21) 

(22) 

(23) 



simplifies the equation for the excitation's amplitude to be found on the n— th edge 

0, (24) 



which is the discrete nonlinear Schrodinger equation (DNLS) [20]. An approximate for 
J/v ^ 1 stationary (consistent with the adiabatic approximation) soliton-like solution 
centered at n = no has the squared amplitude 



V 

8J^ 



-sech 



v{t^ - no) 
4J 



(25) 



Here is the probability density to find the quantum of vibrational excitation on disk 
n. Notice that the evolutionary equation which obeys is the DNLS which has both 
continuous and discrete spectra of solutions. The solution (25) belongs to the discrete 
spectrum of (24) and has a bell-like shape with exponentially vanishing tails. The most 
probable place to find the excitation is in a region around no whose extent is determined 
by v/AJ, that is by the physical parameters in the model: the distance d between the 
disks, the strain tensor constant a, the coupling constants x<f> and xe, and the hopping 
amplitude J. 

Using (21) and imposing the physically consistent conditions rf^x^/a ^ 1 and 
d'^Xe/c( ^ 1 we have the estimate 

< |0n+l - 0n| < 1, < " ^n| < 1- (26) 

Now we turn to the corresponding conformation of F. The total curvature at each disk 



IS 



1 
d 



n\2 



(27) 



sin^ 9- 



nlS^n+l 



+ (^n+1 - e^y 
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and the torsion having fixed the gauge freedom is 

Tn = ""^^"^'^y^"'^"^ (28) 



Taking into account ( 26 ) we may approximate cos ~ 1 and sin 6'n ~ to arrive at an 



approximate expression for the conformation of V in terms of the excitation's presence: 

(29) 



a 



'^dx4> I ,2 
kill 



a 



I 1 2 k 

where |gn| is given by (25). The physical parameters in the model J,a,xeiX<t) 



d establish proportionality between the material response of the a— helix in terms of 
the curvature and the torsion at disk n and the probability amplitude of 
the localization of the quantum of vibrational excitation on that same disk. This 
localization leads to an elastic response in the corresponding configuration of F, with the 



respective energy given by ( 10 ). Provided the five physical parameters are estimated, the 
conformon's energy can be determined and a proper experiment to detect it proposed. 
This is beyond the focus of the present study, although we have a reasonable guess for 
the numerical values of the parameters. Since d is the distance between two adjacent 
disks in the model, it is of the order of the period of the a— helix: d ~ 0.5 nm. The strain 
tensor constant a takes into account the simultaneous response of three hydrogen bonds. 
The spring constant k for the hydrogen bond in the a— helical sections of proteins is 
estimated to be A; = 19.5 N/m [8], leading to an estimate of a ~ 1 x 10~^^ eV.m^|§] 
The hopping amplitude for the Amid-I quantum of vibration along the a— helix is of the 
order of J ~ 1 x 10~^eV [3], [16]. The constants xe and x<t> coupling the conformation of 
r to the on-disk energy are small, most likely smaller (rotational stiffness of chemical 
bonds is smaller that compressional stiffness) that the coupling between phonons and 
the Amid-I quanta (1 x 10"^ eV in the case of the Davydov soliton [21]) thus rendering 



the estimate (26) physically correct. 



4. Conclusions 



We have showed that the vibrational energy localized on the helix acts through a 
curvature and torsion coupling effect to distort the structure of the helix, and then 
the helical distortion reacts again through the curvature and torsion couplings to trap 
the oscillation energy and prevent dispersion. This is the self-localization or self- 
trapping mechanism behind the conformon's existence. Even more surprising is the 
observation that the curvature and torsion of the discrete space curve F used to model the 

§ This estimate can be obtained by equating two equivalent expressions: on one side, the elastic energy 
written for 3 equivalent springs for 10% of the inter-disk distance d deviation from their equilibrium 
position (a non-melting condition according to the empirical Lindemann criterion '22J) and, on the 
other side, the elastic energy written in terms of strain tensor formalism. 
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a— helix are related to the squared probability amplitude of the quantum of vibrational 
excitation, as expressed (29). This shows that the conformation of the a— helix itself is 
quantized, justifying the conformon as a quantum of conformation. 

In this paper we have investigated the role of excitation-conformation interaction on 
the formation of solitary wave packets of conformational energy in a— helical proteins. 
The results obtained here are valid at zero temperature. However, the physical meaning 
of (25) and (29) is experimentally relevant. Pump-probe experiments revealed that 
low frequency nonlinear modes are essential for functionally important conformational 
transitions in proteins containing a— helices [23]. The characteristic lifetime of these 
states is 15 ps (although for the vision-relevant bacteriorhodopsin it can go over 500 
ps). The same experiments also revealed that amino acids and predominantly /3— sheet 
proteins do not have such long-lived states. The proposed model already makes the 
qualitative distinction between sections of the secondary structure of proteins due to the 
full use of the geometry of the hydrogen bonds' network in a— helices. Quantitatively the 
model's viability would be decided upon determination of the physical constants which 
remains an open problem. Other nonlinear mechanisms of energy transport in a— helical 
proteins fail to exhibit such a long lifetime [HI IS] • Thus the conformon proposed here 
might have the answer for the mechanism of direct coherent flow of conformational 
energy for a variety of vital biological processes ranging from electron transfer [21] to 
enzyme action [25] . 
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